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We analyze a new kinematic survey that includes accurate proper motions derived from SDSS DR7 positions, combined with multi- 
epoch measurements from the GSC-II database. By means of the SDSS spectro-photometric data (effective temperature, surface 
gravity, metallicity, and radial velocities), we estimate photometric parallaxes for a sample of 27 000 FGK (sub)dwarfs with [Fe/H]< 
-0.5, which we adopted as tracers of the seven-dimensional space distribution (kinematic phase distribution plus chemical abundance) 
of the thick disk and inner halo within a few kiloparsecs of the Sun. 

We find evidence of a kinematics-metallicity correlation, 3(V (S )/(9[Fe/H]« 40 -H 50 km s -1 dex~', amongst thick disk stars located 
between one and three kiloparsecs from the plane and with abundance -1 <[Fe/H]< -0.5, while no significant correlation is present 
for [Fe/H]^ -0.5. In addition, we estimate a shallow vertical rotation velocity gradient, 3(V<j}/<9|z| = -19 ± 2 km s kpc~', for the 
thick disk between 1 kpc < |z| < 3 kpc, and a low prograde rotation, 37 ± 3 km s - ' for the inner halo up to 4 kpc. 
Finally, we briefly discuss the implications of these findings for the thick disk formation scenarios in the context of CDM hierarchical 
galaxy formation mechanisms and of secular evolutionary processes in galactic disks. 
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The existence of a thic k disk in our Galaxy was revealed by 
Gilm ore & Reid I (119831) . who analyzed starcounts towards the 
South Galactic Pole. Thanks to the many studies carried out 
since then, the main spatial, kinematic, and chemical features 
of this population are well established. Thick disks have been 
also observed in many disk galaxies dYoachim & Dalcanton I 
2006), and they represent the frozen relics of the first phase s 
' of disk galaxy formation (Freeman & Blan d-Hawthorn 1 12001 . 
However, in spite of the many scenarios proposed until now, the 
origin of this component is still unclear. 

In the context of CDM hierarchical galaxy formation mod- 
\ els, it is possible that thick disks are formed by the heat- 
ing of a pre-existing thin disk through a minor merger (e.g. 
IVillalobos & Helmi"ll2008l). b y accretion of stars from disrupted 
satellites (lAbadi et al. 1 12003). or by the s tars formed in situ from 
gas-rich chaotic mergers at high redshift dBrook et al. 120051) . On 
the other hand, simulations suggest that thick disks could simply 
be produced thr ough secular radial migration of stars induced b y 
the spiral arms (Ros kar et al. l2008tlSchonrich & Binnevl2 009). 

In any event, most astronomers agree that our thick disk is 
form ed of an old st ellar population with an age of 8-12 Gyr 
(e.g. lHavwood Il2008l and references therein). The bulk of the 
thick disk stars have metallicity in the range -1 $[ Fe/Hj^ -0.3 ( 
[Fe/Hl^ -0.6, on average ) with enhanced [a/Fe] dBensbv et al. I 
l2005tlReddvetal. Il2006h but note that tails w ith metal-poor 
stars down to [Fe/H]^ -2 dChiba & Beers Il2000t) and metal-rich 
stars up to [Fe/H]^ dBensbv et al. 1120071) have a lso been re- 
vealed. Moreover, according to llvezic et al. I (120081) . a mild ver- 
tical metallicity gradient shifts the mean metallicity to [Fe/H]^ 
-0.8 beyond |z| > 3 kpc. 



The spatial distribution is usually modeled with a symmetric 
exponential density distribution as a function of galactocentric 
coordinates (R,z). Its scale height spans a wide range of mea- 
surements, between h z = 640 pc and 1500 pc, while the local 
normalization varies be etween 13% and 2% in anticorrelation 
with h z (see Fig. 3 of Arnadottir et al. 2008). The distribution 
above the galactic plane is supported by a vertical velocity dis- 
persion, cr w 40 km s _1 , which is associated with an asymmet- 
ric drift of ~ 50 km s , relative to the local standard of rest. 

Significant asymmetries have also bee n detected, such as 
the prominent H ercules thick disk cloud (Par ker et al. I [2003: 
lJuric et ai~ll2008l) . which could correspond to a me rger remnant 
or indicate a triaxial thick disk dLarsen et al. 1120081) . 

In this letter, we present new results regarding the vertical ro- 
tation gradient and, for the first time to our knowledge, evidence 
of a metallicity-rotation correlation in the thick disk stellar pop- 
ulation. 



2. The SDSS - GSC-II catalog 

This study is based on a new kinematic catalog derived by as- 
sembling the astrometric parameters extracted from the database 
used for t he construction of the Second Guide Star Catalog 
(GSC-II; lLasker et al. 1 120081) with spectro-photometric data 
from the Seventh D ata Release of the Sloan Digital Sky Surve y 
(SDSS DR7; e.g. lAbazaiian et all 120091; lYannv et al. 1 12009). 
The SDSS-GSC-II catalog contains positions, proper motions, 
classification, and ugriz photometry for 77 million sources down 
to r ~ 20, over 9000 square-degrees. 

Proper motions are computed by combining multi-epoch po- 
sitions from SDSS DR7 and the GSC-II database. Typically, 5- 
10 observations are available for each source, spanning ~ 50 
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years. Total errors are in the range 2-3 mas yr for 16 < 
r < 18.5, comparable with those of the SDSS proper mo- 
tions dMunn et al. 1 12004). as confirmed by external comparisons 
against QSOs. The construction and properties of this catalog are 
described in detail by Smart et al. (201 0, in preparation), while 
a concise description can be found in Spagn a et al. I J2009). 

Radial velocities {cr Vr < 10 km s -1 ) and astrophysical pa- 
rameters (cr Teff - 150 K, cr log g - 0.25, <t [Fe/H] - 0.20) are avail- 
able for 151000 sources cross-matched with the SDSS spec- 
troscopic catalog. From this list, we select sources with 4500 
K< r eff < 7500 K and logg > 3.5, corresp onding to FGK 
(sub)d warfs, and apply the color thresholds from lKlement et ai~l 
(2009) in order to remove turn-off stars. 

Spectro-photometric distances are computed by means of 
metallicity-depe ndent absolut e magnitude relations, M r - f(g - 
i, [Fe/H]), from llvezic et al. I d2008l) . Here, the observed mag- 
nitudes are co rrected for inte r stellar absorption via the extinc- 
tion maps of Schlegel et al. I d 1998b . while the spectroscopic 
[Fe/H] is used, inst ead of the photometric metallicity applied by 
llvezic et al. I(l2008l) . 

The mean distance of the sample is ~ 2 kpc, while most 
(92%) of the sources are distributed between 0.5 kpc < |z| < 3.5 
kpc and 6 kpc < R < 1 1 kpc. The typical accuracy of the M r cal- 
ibration is 0.3 mag (random) and 0.1 mag (systematic), which 
corresponds to distance errors of Ad/d - 15% and 5%, re- 
spectively. Finally, 3D velocities in the galactocentric reference 
frame, (Vr, V<p, V z ), are de rived by assuming _ 8 kpc, solar 
motion (U Q , V Q , W B ) from lDehnen & Binnev I (G~998), and local 
standard of rest velocity of 220 km s _I . 

In order to produce an accurate sample, we select only stars 
with (f) proper motion errors < 10 mas yr per component, (if) 
errors on the velocity components < 50 km s _1 , (Hi) total ve- 
locity < 600 km s _1 , (fv) distance < 5 kpc, and (v) magnitude 
13.5 < g < 20.5. Overall, the kinematic catalog contains 46000 
stars; in the following sections a subsample of 27 000 low metal- 
licity dwarfs with -3 <[Fe/H]< -0.5 will be used as tracers of 
the inner halo and thick disk and analyzed in details. 

3. Analysis and results 

3.1. Vertical rotation gradient 

Figure Q] shows the distribution of 6538 stars with 1.0 kpc 
< |z| < 1.5 kpc and [Fe/H]< -0.5. In this sample, the contamina- 
tion of thin disk stars is expected to be negligibl^E so that we fit 
the distribution with only two gaussian populations, correspond- 
ing to the thick disk and halo. The least-squares solution of the 
two-component model is good, although the counts at ~ 220 
km s are slightly underestimated (~ -16%) and the velocity 
peak is overestimated of about 7%; this explains a non-optimal 
xl = 3.18. (If we force a third gaussian component correspond- 
ing to the thin disk, the formal goodness of fit improves signifi- 
cantly, xl - 1-37, but the solution becomes ill-conditioned with 
an inaccurate thin disk normalization of (19 + 6)%. ) 

The same procedure is repeated for six height bins: A \z\ = 
0.5-1.0 kpc, 1.0-1.5 kpc, 1.5-2.0 kpc, 2.0-2.5 kpc, 2.5-3.0 kpc, 

1 Assuming a standard model with a thin disk and a thick disk having 
scale-heights of 300 pc and 900 pc, respectively, and a thick disk nor- 
malization of 10% at z=0 pc, about half of the stars belongs to the thin 
disk for 1.0 kpc < |z| < 1.5 kpc, but only a few percent of them with 
[Fe/H]< -0.5 (cfr. e.g.[A umer & Binnev 1 12009). Also, we estimate the 
contamination of metal poor thin disk stars does not exceed 10%, even 
if we adopt a th ick disk with a sh orter h- = 580 pc and a local normal- 
ization of 13% dChen et al. 120011) . 
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Fig. 1. Histogram of the velocity distribution, V$, of the kinematic sam- 
ple with [Fe/H]< -0.5, between |z| =1.0 kpc and 1 .5 kpc. The thick solid 
line shows the best fit of a two Gaussian component model, representing 
the thick disk and halo populations (thin lines). 

Table 1. Parameters of a two-component Gaussian best fit (thick 
disk and halo) for six height intervals. 



Thick Disk Halo 



<|z|> 


N 


<vv> 




<v # > 




/'TP 

ptot 


xl 


(kpc) 




(km s 




(km s 


-') 


(%) 




0.76 


7022 


186+1 


34+1 


46+5 


92+3 


74±2 


5.18 


1.24 


6538 


173+1 


39+1 


32±4 


90+3 


68±2 


3.18 


1.73 


4753 


163+1 


44±2 


29+11 


96+5 


60±4 


1.42 


2.23 


3044 


155+2 


47+3 


36+10 


90+4 


48±5 


1.83 


2.73 


1637 


144+4 


42±5 


49+6 


97±3 


29±6 


1.14 


3.36 


988 


166=fcl 1 


41 + 11 


44+8 


90+4 


13+7 


0.80 



and 3.0^+.0 kpc. The results are reported in TableQ] which lists 
mean height, number of stars, mean rotation velocities and dis- 
persions, fraction of thick disk stars, and reduced^. The halo 
parameters appear quite stable: on average, V$ — Yl ±3 km s _1 
(1 < |z| < 4 kpc), which indicates a slow pro g rade rotation of 
the inner halo, in agreeme nt with some authors (IChiba & Beers I 
120001: iKeplev et all 120071) but different from oth e rs that favor 
a non-rotating inner ha lo (IVallenari et al. 1 120061 ISmith et al. I 
2009; iBond et aL~l 12009). The halo velocity dispersion also ap- 
pears rather constant up to \z\ - 4 kpc, with a mean value of 
cr Vt f, =93+2 km s _1 (uncorrected for the velocity errors). 
Conversely, the thick disk shows a monotonic decreasing of the 
rotation velocity from = 186 kms -1 to 146 km s _1 , for height 
from 0.5 kpc to 3 kpc. In the highest bin (3 kpc < \z\ < 4 kpc), 

increases to 166+11 km s , but we think this is a spurious 
effect of both the larger velocity errors and the small fraction, 
(13+7)%, of thick disk stars that are strongly entangled with the 
halo population. Similarly, in the same z-range, the velocity dis- 
persion increases from cry^ = 34 km s _1 to ~45 km s _1 , in part 
because of the tangential velocity errors that scale with distance. 

We exclude the highest bin and also the lowest, as it is prob- 
ably contaminated by thin disk stars which are difficult to de- 
convolve from the thick disk population. Thus, we estimate the 
gradient, 

d(V^)/d \z\ = -19 + 2 km s" 1 kpc" 1 (1) 
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Fig. 2. Velocity-metallicity distribution of 20251 stars with |z| =1.0-3.0 kpc and [Fe/H]< -0.3. The dashed line indicates the thick 
disk rotation, = 173 km s at (|z|) = 1.24 kpc (Table [TJ. The box defines the region, shown in Fig. [3] in which the thick disk 
population dominates. 




Fig. 3. Iso-density contours of 13 108 relatively metal-poor stars with —1.0 <[Fe/H]< -0.3 and \z\ =1.0-3.0 kpc. White crosses mark the ridge 
line of the maximum likelihood V$ vs. [Fe/H]. 
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Fig. 4. vs. [Fe/H] distribution of stars with -1.0 <[Fe/H]< -0.5. Black crosses mark the rejected stars beyond 3cr of the velocity 



ellipsoid of the thick disk to minimize the contamination from halo stars (see Sect. l3.21 >. The solid line connects the mean velocities 
(V^), which are plotted with 2<x error bars. The dashed lines indicate the +lcr spread of the velocity distribution. Top, middle, and 
bottom panels refer to |z| =1.0-1.5 kpc, 1.5-2.0 kpc, and 2.0-3.0 kpc, respectively. 
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Table 2. Kinematics-metallicity correlation of thick disk stars 
with -1.0 <[Fe/H]< -0.5 for three height intervals. 



<|z|> 


N w 




<9<^>/d[Fe/H] 


Ps 


(kpc) 






(km s -1 dex~') 


(xlO- 2 ) 






3cr 2cr 


3tr 2cr 


3cr 2cr 


1.23 


3994 


3672 2915 


50±5 39+5 


17±2 15±2 


1.73 


2641 


2348 1715 


54+6 35+5 


18±2 16±2 


2.37 


2194 


1768 1131 


35±8 33±5 


10±2 14±3 



and the extrapolated intercept, V^iz = 0) = 196 + 3 km s _1 . 
Our result is significantly smaller than the value , -30 ± 3 
km s _1 kpc -1 , measured by Chiba & Beers I d2000l) . who ana- 
lyzed stars with abundance in the range, -0.8 <[Fe/H]< -0.6, 
where the thi c k disk dominates. A similar tr end was estimated by 
iGirard et al.l d2006l) . ICarollo et all (120091) . and by iBond et al. I 
(2009), although they adopted a nonlinear function. 

Instead, a shallower slope was found by iMaiewskil (119921) . 
who derived a gradient of -21 + 1 km s -1 kpc -1 for |z| < 5 kpc, 
after separating the halo population from tha t of the thick disk. A 
low kinematical gradient was also found by Spag na et al. I {l996) 
and, more recently, by All ende Prieto et al. I (120061) . who esti- 
mated -10 km s kpc 1 and -16 km s _1 kpc -1 , respectively. 
The difference between these results can be explained, at least in 
part, by thin disk and halo stars contamination, which tends to 
produce steeper velocity gradients. 



3.2. Rotation - metallicity correlation 

The disk and halo populations are apparent in the vs. [Fe/H] 
distribution (see Fig. [2j». In particular, the region -1.0 <[Fe/H]< 
-0.5 and km s _1 < V$ < 300 km s _I does contain the bulk 
of the regular thick disk stars, besides a small number of stars 
belonging to the metal-poor tail of the thin disk and to the high- 
metallicity tail of the inner halo. Actually, a significant fraction 
of thin disk stars are expected for [Fe/H]> -0.5, while towards 
lower abundances, [Fe/H]< -1, the thick disk metal weak tail 
and th e newly discovered flattened inner halo (Mor rison et al. I 
2009) are also present. 

Figure [3] shows the iso-density contours of the velocity- 
metallicity distribution of stars with \z \ = 1.0- 3.0 kpc and 
-1.0 <[Fe/H]< -0.3. As in llvezic et all J2OO8I) and lBond et al. I 
(120091) . no correlation appears in the transition region between 
the thin and thick disks ( [Fe/H]S -0.5). Instead, we notice a 
shallow but clear slope for [Fe/H]^ -0.5, undetected by pre- 
vious studies, which indicates that the metal-rich stars tend to 
rotate faster than the metal-poor ones. In particular, the top- 
density ridge increases from V$ - 150 km s _I at [Fe/H] « -1 
to V = 170 km s" 1 at [Fe/H]*: -0.4. Inspection of Fig.[3]also 
proves a bimodal distribution with a secondary maximum lo- 
cated at [Fe/H]s; -0.55, close to the value of the mean metallic- 
ity of the thick disk, and the peak at [Fe/H]^ -0.38 due to thin 
disk stars. 

To quantify the correlation, we first select the stars within 
A[Fe/H]=0.05 bins in the range -1.0 <[Fe/H]< -0.5 and located 
at the different height intervals: A |z| = 1.0-1.5 kpc, 1.5-2.0 kpc, 
and 2.0-3.0 kpc. Then, the stars with velocities (Vr, V^, V z ) out- 
side 3cr from the thick disk velocity ellipsoid, corresponding to 
a confidence level of 97.1%, were rejected to minimize the con- 
tamination from the halo stars. We adopted (V^) as a function 
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Fig. 5. V distribution for [Fe/H]< -0.5, -0.7, and -1.0, in de- 
scending order. The solid lines show the best fit of two Gaussian- 
component models (thick disk and halo). Bottom, middle, and 
top panels refer to \z\ =1.0-1.5 kpc, 1.5-2.0 kpc, and 2.0-3.0 
kpc, respectively. The solid and dashed lines mark the values 
V 6 = km s _1 and 170 km s _1 . 



of z derived from Table Q] and assumed constant dispersions: 
o- v< p = crvz = 40 km s _1 and ct Vr = 1.5 cry^ = 60 km s _1 . 

Finally, mean velocities were computed for the bona fide 
thick disk stars and the slope, <9(V0)/3[Fe/H], is estimated by 
means of a linear fit for the height intervals A \z\ =1.0-1.5 kpc, 
1.5-2.0 kpc, and 2.0-3.0 kpc. For each bin, mean height, to- 
tal number of stars, number of stars used (after 3<r and 2<x re- 
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Table 3. Fitted parameters, as in Table[TJ for different metallicity intervals, -3.0 < [Fe/H] < [Fe/H] max , where -1.0 < [Fe/H] max < 
-0.5. 



Thick Disk Halo 



[Fe/H] raax N <V > try, <V„> ov, ^ X \ 

(dex) (kms- 1 ) (kms- 1 ) (%) 



1.0 kpc < \z\ < 1.5 kpc 


-0.5 


6537 


173 


± 1 


39 ± 1 


33 ±4 


90 


±3 


68 


± 2 


3.19 


-0.6 


5470 


170 


± 1 


39 ± 1 


35+7 


94 


± 4 


61 


± 3 


2.13 


-0.7 


4511 


167 


± 1 


39 ±2 


33+7 


94 


± 4 


55 


± 3 


1.80 


-0.8 


3675 


165 


± 1 


38 ± 3 


31 ±6 


93 


± 3 


46 


± 3 


1.57 


-0.9 


3036 


162 


± 2 


38 ±2 


29 ±6 


93 


±3 


38 


±4 


1.51 


-1.0 


2543 


162 


± 2 


37 ±2 


24 ±5 


92 ± 3 


30 


±4 


1.59 


1.5 kpc < |z| < 2.0 kpc 


-0.5 


4753 


163 


± 1 


44 ± 2 


29 ± 10 


96 


±5 


60 


±4 


1.42 


-0.6 


4113 


157 


± 2 


45 ± 2 


23 ± 10 


94 


± 5 


56 


±4 


1.49 


-0.7 


3480 


154 


± 2 


45 ± 2 


23 ± 11 


95 


± 5 


49 


+ 5 


1.28 


-0.8 


2936 


152 


± 2 


43 ± 3 


22 + 9 


95 


± 4 


41 


± 5 


1.30 


-0.9 


2488 


149 


± 3 


42 ± 3 


19 + 8 


94 


± 4 


33 


± 5 


1.20 


-1.0 


2112 


145 


±4 


44 ±5 


18+9 


94 


±4 


25 


±6 


1.18 


2.0 kpc < \z\ < 3.0 kpc 


-0.5 


4680 


152 


± 2 


45 ± 2 


43 ±6 


94 


± 2 


40 


±4 


1.72 


-0.6 


4176 


150 


± 3 


46 ±3 


42 ± 6 


94 


± 2 


34 


± 4 


1.56 


-0.7 


3671 


147 


± 3 


45 ± 3 


37 ±6 


93 


±2 


30 


±4 


1.41 


-0.8 


3183 


144 


±4 


43 ± 4 


33 ±5 


92 


± 2 


25 


± 4 


1.31 


-0.9 


2830 


143 


±4 


40 ± 5 


32 ±5 


92 


± 2 


19 


± 4 


1.20 


-1.0 


2486 


141 


±5 


38 ± 6 


30 ±4 


93 


±2 


14 


±4 


1.20 



jection), slope, and Spearman's rank correlation coefficient are 
listed in Table f2] while the observed distributions are shown in 
Fig. |4] Overall, a kinematic-metallicity correlation of about 50 
km s" 1 dex 1 is detected up to \z\ - 2 kpc, while a shallower 
slope (~ 35 km s _1 dex -1 ) is present between 2 < |z| < 3 kpc. It 
is possible that these values are affected by a residual contamina- 
tion of halo stars, whose presence can be inferred by the number 
of rejected high velocity stars shown in Table [2] being greater 
than the 3% expected in the case of a pure Gaussian distribu- 
tion. Nevertheless, even if we apply a conservative 2<x selection 
(73.8% confidence level), we still find a correlation at the level 
of 30H-40 km s dex -1 , as reported in the last column of Table 

This conclusion is consistent with the systematic slowing 
down of the thick disk rotation, which results from fitting a 
two Gaussian-component model, representing the thick disk and 
halo populations, as more metal-poor thresholds are applied: 
[Fe/H] max < -0.5, < -0.6, ... < -1.0 (see Table©. This effect 
is depicted in Fig. [5] which shows how the thick disk compo- 
nent both decreases and shifts towards lower V# values, when 
different subsamples of metal poor stars are selected. 

In addition, we estimate the rotation-metallicity correlation 
by fitting the thick disk ( V^) values from Table [3] through the 
following integral linear model: 

<V*> = We/H] ) + a ■ «[Fe/H]> - [Fe/H]„) , (2) 

where ([Fe/H]) is the average for the stars with -3 < 
[Fe/H]<[Fe/H] max , a = d(V+)/d[Fe/H\, and V,([Fe/H] ) is the 
mean velocity of the reference metallicity, which we set to 
[Fe/H]o = -0.6. In Figure [6] the lines connect the values from 
Eq.|2]at the different [Fe/H] max thresholds. These results confirm 
both a vertical gradient consistent with the value derived in Sect. 
13.11 and a rotation-metallicity correlation in the range of 40H-50 
km s _1 dex 1 for the thick disk. 



We also considered the hypothesis that a false trend 
vs. [Fe/H] might derive from the tangential velocity esti- 
mated through the metallicity-dependent photometric parallaxes. 
Actually, the correlation would still be significant even if the 
M,-calibration were subjected to a systematic error up to 0.4 
mag per dex. Moreover, no kinematics-metallicity correlation is 
expected to arise because of the color-selection criteria of the 
SDSS spectroscopic targets, which although they produce a bias 
towards metal poor stars, cannot affect the conditional prob- 
ability distribution at a given metallicity, Pr(V0 I [Fe/H]), and no 
further kinematical selection is applied. 

Thus, we conclude that the observed correlation is an intrinsic 
signature of our sample. 

4. Discussion and conclusions 

The existence of a vertical velocity gradient and a rotation- 
metallicity correlation sets important constraints on the origin of 
the thick disk. The estimated gradient of -19 ± 2 km s kpc -1 
is consistent with Nbody simulations of disks thickened by a 
singl e minor merger with a low/intermediate orbital inclination 
(e.g. IVillalobos & Helrriil 120081) . as well as by the interaction 
with n umer ous dark subhalos, a s discu ssed bv lHavashi & C hiba 
( 2006j) and iKazantzidis et al. I d2008l) . whose simulations show 
kinematic gradients of -(10 4- 30) km s _1 kpc -1 and of -20 
km s _I kpc -1 , respectively, for 1 kpc < \z\ < 3 kpc. A vertical 
rotation gradie nt of about -20 km s -1 kpc can also be inferred 
from Fig. 5 of lAbadi et al. I (|2003), who investigated thick disks 
formed by accretion of both the stars of a pre-existing thin disk 
and the debris from disrupted satellites. Unfortunately, we have 
not found any explicit kinematic predic tion in the scenario of the 
chaotic gas-rich mergers described by [Brook et al. I (120051) . al- 
though lHavashi & Chiba I (120061) state that a velocity shear "may 
have difficulties in this regard". Finally, to the best of our knowl- 
edge, explicit predictions of kinematics-metallicity correlations 
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Fig. 6. Mean vs. [Fe/H] max for different abundance ranges, 
-3 <[Fe/H]<[Fe/H] milx , with -1.0 <[Fe/H] raax < -0.5. The three lines 
show the (V^> computed from Eq.[2]with fitted slopes, 3<V^>/3[Fe/H] = 
(39 ± 5, 56 ± 9, 38 ± 13) km s _1 dex -1 , and zero-points, V>([Fe/H] ) = 
(178±1, 170±2, 158±3)kms-» for \z\ =(1.0-1.5, 1.5-2.0, 2.0-3.0) kpc, 
respectively. 



are missing in the current CDM scenarios of satellite accretion 
or minor mergers. Hopefully, our results will motivate theoreti- 
cians to investigate this issue in their future models. 

In the context of models based on disk secular pro- 
cesses of stellar migration driven by interactions with spi- 
ral arms, a vertical gradient of ~ -15 km s _I kpc 1 is re- 
ported by Loebman et al. (2008), who, conversely, did not de- 
tected any vs. TFe/Hl correlati on. The simulations carried 
out bv lSchonrich & Binnevl d2009l) indicate a mild trend (~ 10 
km s _1 dex~') at z ~ kpc, which decreases with height and dis- 
appears for |z| ^ 1 kpc. Possibly, by adopting appropriate param- 
eters, their inside-out disk formation model could reproduce the 
observed downtrend (Schonrich, 2009, private communication). 
Thus, more attention should be devoted to this scenario as a pos- 
sible theoretical framework to explain the rotation-metallicity 
relation in the thick disk of the Milky Way. 
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